Instabilities of the AA-stacked graphene bilayer 
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Tight-binding calculations predict that the AA-stacked graphene bilayer has one electron and one 
hole conducting bands, and that the Fermi surfaces of these bands coincide. We demonstrate that 
as a result of this degeneracy, the bilayer becomes unstable with respect to a set of spontaneous 
symmetry violations. Which of the symmetries is broken depends on the microscopic details of 
the system. We find that antiferromagnetism is the more stable order parameter. This order is 
stabilized by the strong on-site Coulomb repulsion. For an on-site repulsion energy typical for 
graphene systems, the antiferromagnetic gap can exist up to room temperatures. 

PACS numbers: 73.22.Pr, 73.22. Gk, 73.21.Ac 



Introduction. — Graphene is a zero-gap semiconductor 
demonstrating a host of unusual electronic properties [D- 
[jf]. In recent years, the synthesis of bilayer graphene trig- 
gered investigations of the bilayer systems, partly driven 
by the desire to create graphene-based materials with 
an electron gap. Moreover, the graphene bilayers are 
interesting materials in their own right. Most efforts 
have focused on the study of the AB-stacked bilayer 
Q for which high-quality samples are available 0, Q . 
Lately, the experimental realization of the AA-stacked 
graphene has been reported @, HI ■ In this paper we dis- 
cuss electronic properties of the AA-stacked graphene bi- 
layer (AA-BLG), which, until recently, received very lim- 
ited theoretical attention 

It is known that the AA-BLG tight-binding spectrum 
has four bands, of which one electron band and one 
hole band cross the Fermi energy [lOfl . The Fermi sur- 
faces of these two bands coincide [3, This fea- 
ture has drastic consequences for the electronic prop- 
erties of the bilayer because it enables several electron 
and electron-phonon instabilities, including: antiferro- 
magnetism (AFM), current-ordered states, bilayer exci- 
ton condensation, and instability toward the shear shift 
of the layers. The type of ground state order depends on 
the microscopic details of the system and can be changed 
by applying stress, external pressure, the presence or ab- 
sence of the substrate, etc. Below we will limit our at- 
tention to the AFM order and the structural instability 
with respect to the shear layer shift (shear instability 
for short). These two choices are justified. The on-site 
Coulomb repulsion is the strongest interaction in the AA- 
BLG system, and this interaction is sufficient to guaran- 
tee the stability or mctastability of the AFM order. As 
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for the shear instability, there are experimental 
and numerical [15[ suggestions that AA-stacked graphene 
multilayers may be unstable with respect to the mechan- 
ical displacement of the layers with respect to each other. 
However, our calculations show that the shear instabil- 
ity driven by the conducting electrons seems to have a 



crossover temperature which is too low to be experimen- 
tally observable. 

The model. — In the AA-BLG, carbon atoms of the 
upper layer are located on top of the equivalent atoms of 
the bottom layer. The system is modeled by the tight- 
binging Hamiltonian for p z electrons of carbon atoms 
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nihilation operators of an electron with spin a in the layer 
% = 1, 2 on the sublattice A (B) at site n E A (m £ B). 
The amplitude t (to) in Eq. (|TJ) describes the in-plane 
(inter-plane) nearest-neighbor hopping, while t g corre- 
sponds to the inter-layer next-nearest neighbor hopping. 
The interplane distance in bilayer graphene c ~ 3.3 A [8(], 
and it is larger than the in-plane carbon-carbon distance 
a ~ 1.4 A. Thus, the in-plane hopping integral t is larger 
than the interplane one, to. For calculations we will 
use the characteristic values t w 2.57 eV, to ~ 0.36 eV, 
t g « — 0.03eV [l6[. We omit next next-nearest neigh- 
bor hopping between A (B) sites since the corresponding 
term only shifts the zero-energy level. 

The elementary unit cell of bilayer graphene consists 
of four atoms. It is convenient to introduce the bi- 
spinors ip^ = ^V'k*4o-> V'kScr)' w ^h spinor components 

^Ia* = ("La' a La) and ^Ib* 

where <^ k = arg{/ k }, and 
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FIG. 1: (Color online) (a) The band structure of the AA- 
stacked bilayer graphene. (b) The k-dependence of the spec- 
tra e^ 1 near the Dirac point K, located at momentum K; 
k = K + Sk y e y . Bands s — 2 and s = 3 intersect at the Fermi 
level £f. (c) Solid (green) lines show six arcs of the Fermi 
surface in the first Brillouin zone. 



Let us define a set of Pauli matrices f a acting on the 
layer index, and a second set a a of Pauli matrices acting 
on the sublattice index. In terms of these matrices, the 
bilayer Hamiltonian in k-space can be written as i?ok = 
- [toT x + (t + t g T x )a x \f k \], or, explicitly 
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The Hamiltonian ([3]) is invariant under the transposi- 
tion of the sublattices and of the graphene layers. That 
is, [a x , if k] = [t x , ffok] = 0. Thus, the eigenvectors of 
the matrix ([3]) can be classified according to the quantum 
numbers a and t, which characterize the eigenvector par- 
ity under a x and t x transformations. Using these sym- 
metries it is easy to find the transformation which diago- 
nalizes H^: it is U = (r x + f z ) (a x + a z ) /2 = ?7 _1 . The 

(s) 

electron spectrum obtained consists of four bands, 
and each band has a unique value of the pair {a, r): 
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(7 = 1, T = 1, (4) 

a=l, r=-l, (5) 
a = -l, r=l, (6) 
<r = -l, r = -l. (7) 



The band structure is shown in Fig. [T] The bands s = 2 
and s = 3 cross the Fermi energy level near the Dirac 
point K,, located at momentum K = 2w{y/3, l}/(3\/3a) 
and the Dirac point K! located at momentum K' = 



2tt{\/3, -l}/(3 v / 3a) [see Fig. [H»]. The most inter- 
esting feature of this band structure is that at half fill- 
ing (which corresponds to undoped AA-BLG) the Fermi 
surfaces of both bands coincide. The Fermi level is 
£f = tgto/t, while the Fermi surfaces are given by the 
equation |/k| = to/t. For to/t -C 1 one can expand the 
function |/k| near the Dirac points and demonstrate that 
the Fermi surface consists of six arcs inside the first Bril- 
louin zone with the radius k r = 2£o/(3ta) [Fig. He)]. 

The matching of the Fermi surfaces turns out to be 
quite stable against changes in the tight-binding Hamil- 
tonian. First, it survives if we add more distant hopping 
terms to Hq. Moreover, even layer-asymmetric systems 
(e.g ., similar to the single-side hydrogenated graphene 
17. 1) may posses this property. However, it is clear that 
the different types of interactions, e.g., electron-electron 
or electron-phonon ones, can destabilize such a degener- 
ate spectrum. 

Mean-field Hamiltonian. — The presence of two bands 
with identical Fermi surfaces makes the system unstable 
with respect to spontaneous symmetry breaking. We will 
demonstrate that the Hamiltonian symmetries a x and t x 
can be used to narrow the possible symmetry choices. 

In the mean-field approach, the two-particle interac- 
tion operator H[ nt oc xp'tp'tpip is replaced by a single- 
particle operator SH lnt oc (iflijjjil^ij), where the aver- 
age (ijj^ijj) is different types of non-superconducting or- 
der parameter. The values of these order parameters are 
found from the self-consistency conditions. To be at least 
metastable, the order parameter must open a gap at the 
Fermi level. The most general form of 5Hi nt , which can 
open an insulating gap, is 



SH, 



(8) 



SH k 



where A°i , 
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j , and A£f are real- valued order pa- 
rameters, which, in general, are functions of k. To open 
a gap, the corresponding term in SHi nt must couple the 
conducting bands and s£ 3 . Since these bands have 
unequal values of a and r [sec Eq. (J5J) and Eq. , there- 
fore, only terms containing A£f with a, ^ x may cou- 
ple these bands. Other terms commute either with a x or 

Tx- 

To find the renormalized spectrum of the bands near 
the Fermi-level, we should diagonalize the matrix H^ = 
Hok + SH^. Performing the unitary transformation of 
-ffker with U written above, we obtain 
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Other elements of this matrix arc unimportant for further 
consideration. Solving the secular equation for the 2x2 
matrix presented in Eq. we obtain the renormalized 
spectrum of the bands with s = 2,3: 
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The gap between renormalized bands is equal to Ao = 
2rriin k |A kcr |. We see that the contribution to the gap 



comes only from A k 2 , Aj^, A^., 



and A k order param- 



eters which break down both sublattice and layer sym- 
metries. It is easy to show that other elements of the full 
4x4 matrix © give only the second-order contribution 
to this result. Taking other order parameters zero, the 
matrix SH^ can be written in the form 
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In principle there is a huge number of possible insta- 
bilities generated by different types of interactions which 
can lead to the mean-field interaction Hamiltonian of the 
form in Eq. (JTTJ) . Below we will show that A k 2 can be 
related to the G-type antiferromagnetic order parameter 
(i.e., each spin is antiparallel to all nearest-neighboring 
spins) produced by the on-site Coulomb repulsion. The 
A k ^ can be attributed to the instability toward the ho- 
mogeneous shift of the one graphenc layer with respect 
to another one. The order parameters Aj^ and can 
correspond, e.g., to excitons, which produce a current 
flowing inside and between the layers, respectively. 

Antiferromagnetic state. — It is known that the 
Coulomb interaction among electrons in graphene is 
rather strong and the value of the on-site Coulomb re- 
pulsion energy U is about 10 eV jla |. However, graphene 
remains semimctal since the electron density of states at 
the Fermi level is zero. In contrast, the AA-BLG has a 
Fermi surface and the density of states at the Fermi level 
is finite. Then, one can expect that the role of electron- 
electron interactions in AA-BLG is more important and 
it can affect the ground state. We restrict ourselves here 
by considering the on-site Coulomb interaction and write 
the Hubbard Hamiltonian in the form 
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where n niAa = a nwT a ni(T , n mlBa 
5" = —a. It is known that the ground state of the 
Hubbard-likc models at half-filling can be antiferromag- 
netic (AFM). For the AA-BLG symmetry, three types 
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FIG. 2: (Color online) The dependence of the AFM gap A 
on the on-site Coulomb repulsion U. Solid (red) curve is cal- 
culated by solving Eq. (|14|l . while the dashed (blue) curve is 
calculated from Eq. (|15[) . The inset shows the electron spec- 
trum near K, point at U = 5.5 eV (A m 0.12 eV). 



of AFM ordering (having different spin arrangement in- 
side the unit cell) are possible. However, only the G-type 
AFM order (AFM arrangement both between sublattices 
and layers) opens a gap at the Fermi level even if the in- 
teraction is arbitrary small. 

In the mean-field approximation we represent n n i Q(T 
(a = A, B) in Eq. (JTSJ in the form n macr = n ia< j + 8n nlatT , 
where n taa = (n ma(T ) and 8n ma<T = n n iac - ri iaa . 
The mean-field Hamiltonian is obtained then by neglect- 
ing the terms quadratic in 8n n i aiJ . For G-typc AFM, 
the spin-up and spin-down electron densities are redis- 
tributed as ni At = n 2 Bf = n 2 Ai = n\Bi = (1 + An)/ 2 
and niAi = n 2Bi = n 2At = n iBf = (1 ~~ An)/2, while the 
total on-site electron density ni aa +ni aS remains equal to 
unity. Thus, the mean-field interaction Hamiltonian has 
the form in Eq. © with A k ^ = —A, A k ^ = +A, where 
A = UAn/2. Other terms in Eq. ([8]) are equal to zero. 



The eigenvalues E^ and eigenvectors fp ko . of the ma- 
trices i? kcr = i?ok + 8Hk a can be found analytically. The 
spectra of spin-up and spin-down electrons are equal. 
The spectra of the bands 2 and 3 have the form in 
Eq. CGI with fe k 2 J = fe k 3 j = 0, and |A kff | = A. The 
spectra of the bands 1 and 4 are given by Eq. ([TO"]) with 
the replacement £ k 2 ' 3 ^ — > £ k • When the gap is open, 
the lower two bands are filled, while the upper two are 
empty. To find the value of the gap, one needs to solve 
the self-consistent equation for An = 2A/U: 
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where Vbz = 8n 2 / (3^/3a 2 ) is the volume of the first Bril- 
louin zone. We introduce the dimcnsionlcss density of 
states po(C) = f dk-8(C — |/k|)/%z, which is related [lj 
to the graphene density of states p gr (E) according to 
Pgr(E) = po(\E/t\)/t. Equation (flU)) then becomes 
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where 5 = A/t and £o = ^oA- The integral of the second 
term in the left-hand side of Eq. (|14[) diverges logarithmi- 
cally when A — > 0. In the limit of small A, from Eq. (TT4")) 
one can derive 



A = 2V*o(3t-t ) 
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Figure [2] shows the dependence of A on U. Taking 
the value of U = 8 -f- 9eV 01, we obtain A = 2 3eV 
and the magnetic moment at each site /xeAn about lpB- 
However, we do not know exact value of U for AA-BLG. 
These calculations were done at zero temperature. At 
finite temperatures no long-range AFM order exists. The 
crossover temperature T* between the short-range AFM 
state and paramagnetic state can be estimated as T* ~ 
A/**. 

Discussion. — Other possible types of ordering could be 
considered following the same approach used for AFM. 
However, whether a particular order is stable and ob- 
servable depends on the values of the hopping amplitudes 
and a characteristic energy of the appropriate interaction. 
For example, applying the mean field approximation to 
the model with the on-site repulsion we find that the 
charge density is unstable for our choice of parameters. 

The next evident possibility to open a gap in the spec- 
trum is to shear one graphene layer with respect to an- 
other. The shift u deforms the shape of the unit cell 
changing AB bonds between different layers, giving rise 
to the appearance of the order parameter A^. Assuming 
that the hopping amplitude t g changes linearly with u, we 
can write for different bonds t g (u) ss t g (0) ± (dt g /du)u. 
Now the electronic energy of the system becomes a func- 
tion of u. Taking into account the elastic contribution 
C s hW 2 /2 (where C s h is the corresponding shear modulus) 
to the total energy and minimizing this energy with re- 
spect to u, we obtain the value of the equilibrium shift 
between the layers 
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If we assume that C s h is approximately equal to the shear 
modulus in graphite and \dt g /du\ ~ t g /a, we conclude 
that the shift u eq and the corresponding energy gain are 
too small to be observable. However, this conclusion 
must be taken cautiously. First, we have no accurate 
information on C s h and \dt g /du\, whose precise values 
are very important for the estimate of u eq and the cor- 
responding energy. Moreover, u cq could be enhanced by 
pressure or the presence of a substrate. Finally, the shift 
can be induced by a different mechanism. 

In conclusion, we demonstrate that the AA-BLG is 
unstable with respect to a set of symmetry-breaking in- 
stabilities, which can give rise to the existence of several 
order parameters of different nature. We show that the 
AFM order can be observed in the system. The pos- 
sible existence of other types of orders in the AA-BLG 
depends on the system parameters and the external con- 
ditions (temperature, pressure, substrate, etc). 
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